3. Results and discussionWe study the thermally activated electrons transport behaviors through the discrete dopant-induced QDs in the impurity band of the silicon JNT within the temperature range from 6 K to 250 K by the temperature-dependent conductance characteristics of the silicon JNT. The thermally activated electron transport is dominated by phonon-assisted hopping through several dopant-induced QDs.[21] Mott believed that a hopping electron with activation energy Ea would always try to find a lower energy state around the Fermi level by the variable range hopping (VRH) in the absence of long-range Coulomb interaction.[22] At high temperatures, the temperature dependence of conductance G for the nearest neighbor hopping (NNH) exhibits as
, where kB is the Boltzmann constant.[23] If the activation energy at low temperatures is reduced as large as the Coulomb interaction energy, the VRH conductance which has been predicted by Efros and Shklovskii (ES) should obey the 1/2 exponent law given by
, where TES is the temperature-independent ES coefficient.[24] Therefore, with temperature increasing, three kinds of electron hopping transports would be observed in the following consequence: ES-type VRH, Mott-type VRH, and NNH.[21] In our experiment, we find two obvious transition platforms within the critical temperature regimes of A and B for the experimental conductance data, which are extracted from the unified transfer characteristics at the gate voltage positions of the initial transconductance gm peak in Vg1 and valley in Vg2. One crossover temperature TA in the higher temperature regime A corresponds to the transition of thermally-activated-electron transport behaviors from NNH to Mott-type VRH.[25] The other crossover temperature TC in the lower temperature regime B corresponds to the VRH transition for the delocalized electrons from Mott law to the ES law under the influence of Coulomb interaction.[26] The quantitative analysis on the experimental temperature-dependent conductance data has been taken by linear fitting of
for NNH,
for Mott-type VRH, and
for ES-type VRH.
Figure 2(a) presents the temperature-dependent conductance G (i.e.,
) curves obtained from the transfer characteristics at low source–drain bias VDS = 1 mV. The clear quantized current steps below the temperature of 75 K indicate the population of the individual sub-bands caused by the quantum confinement effect.[27] Above the temperature of 75 K, the thermal energy is greater than the subbands spacing, resulting in the smearing of the current steps. Figure 2(b) illustrates the corresponding transconductance gm–VG characteristics at different temperatures. We define the onset gate voltage Vgt as the initial point of the first gm peak, where the drain current in the conduction channel is at the onset state. The onset gate voltage Vgt increases from 6 K to 20 K and decreases from 100 K to 250 K, resulting from the interaction of the induced image charges in the dielectric interface with the impurity and the subband states in the channel.[13] To explore the conductance characteristics under the same filled energy level, we take the onset gate voltage Vgt to unify the transfer characteristics in Fig. 2(a) for the alignment of energy levels at different temperatures. In order to clarify the electron hopping transport in the impurity band, figure 2(c) shows the IDS–Vg curves (upper part) and the corresponding transconductance gm–Vg curves (lower part) at the temperature of 6 K under the bias VDS varying from 1 mV to 10 mV. The clear oscillatory current with several splitting peaks (upper part) identifies the coupling of dopant-induced QDs in the impurity band. With the gate voltage increasing, the Fermi energy level of electrons in the quantum confined channel is allowed to enter the conduction subbands, resulting in the current steps.[28, 29]
In order to study the temperature-dependent conductance G characteristics in Fig. 2(a), we firstly extract the experimental conductance data for different temperatures according to the initial gm peak at the gate voltage Vg1, in which the effective mobility of hopping electrons is the highest in the impurity band. For comparison, we also provide the temperature-dependent conductance data according to the initial gm valley at the gate voltage Vg2, in which the electron hopping behavior is suppressed by the Coulomb interaction. Figure 2(d) provides the Arrhenius conductance plots (G vs. 1/T) at the gate voltages of Vg1 and Vg2 within the temperature range from 6 K to 250 K. Interestingly, it is found that the conductance at the temperature of 6 K is much larger than that of 10 K, which may be related with the Coulomb interaction. As the temperature increases, the incomplete ionized donors in the channel are gradually transformed into ionized donors, which result in stronger impurity scattering to reduced the electron mobility.[30] The inset G–T diagrams in Fig. 2(d) show two apparent platforms of the conductance for the gate voltages Vg1 and Vg2 within the temperature regions of around 200 K and around 100 K, as the color marked A and B. The platform of the temperature-dependent conductance for the gate voltage Vg1 within the temperature range from 50 K to 100 K is more evident than that of the gate voltage Vg2 within the temperature range from 75 K to 125 K, resulting from the stronger quantum confinement at the initial stage of the conduction channel in the silicon JNT.
In order to precisely determine the transition temperatures TA in temperature region A, we replot the Arrhenius curves of temperature-dependent conductance at high temperature regime as the curves of ln G vs. 1/T in Fig. 3 for the gate voltages Vg1 and Vg2 respectively. The inset shows the Mott-type VRH temperature dependence of the conductance in ln G vs. (1/T)1/4 scales. Here, we are interested in the hopping behaviors of the thermally activated electrons from the ionized dopant atoms. The linearly fitting by
is shown in Fig. 3 for the thermally activated electrons by NNH above the crossover temperature TA, in which the electrons have enough thermal activation energy Ea to overcome the potential barriers between the nearest neighbor states. The hopping distance r in NNH is supposed to be equivalent to the mean distance d between neighbored dopant atoms, which can be determined by the doping concentration Nd as
According to the doping concentration
, the hopping distance between the nearest neighbored dopant atoms is estimated to be 13.4 nm. According to the linearly fitting of the experimental conductance data in Fig.
3 within the temperature range from 175 K to 250 K, the activation energies
Ea for the hopping electrons are extracted to be 45.2 meV and 33.7 meV respectively for the curves at the gate voltages
Vg1 and
Vg2.
[31] According to Mott’s theory, the activation energy
Ea for electron hopping is related to the constant density of states (DOS)
g0 at the Fermi level and the hopping distance
r as
[32]The probability
P for hopping is proportional to the conductance
G of the device, which depends on the overlap integral of the wavefunctions within the localization radius
a and the activation energy
Ea,
[26] i.e.,
Taking the expression of
Ea (Eq. (
2)) into Eq. (
3), one may get the relation between the probability and the hopping distance,
. After taking the derivative of the hopping probability
, the optimal electron hopping distance
is obtained to be
[33]where
aM is the localization length in the Mott hopping regime. At the crossover temperature
TA, the electron hopping distance of VRH is assumed to be equal to that of NNH,
By taking the expression of
(Eq. (
4)) into Eq. (
5), we deduce the expression
which is determined by the localization length
aM and the density of states
g0. In order to determine the localization length
aM, we may use the Mott-type VRH conductance expression by substituting Eqs. (
2) and (
4) into Eq. (
3), i.e.,
. The temperature factor
TM is related to the localization length
aM by
. According to the fitting results from 175 K to 250 K of ln
G vs.
scales in the inset of Fig.
3, the temperature factor
TM extracted from the slopes is 38.78
4 K for the gate voltage
Vg1 and 28.9
4 K and for the gate voltage
Vg2. Therefore, the localization length
aM can be obtained from the expression of
TM be about 3.47 nm and 4.66 nm respectively for the gate voltages
Vg1 and
Vg2. The density of states
g0 at the Fermi energy can be estimated to be
and
from Eqs. (
2) and (
5) by taking the activation energies of 45.2 meV and 33.7 meV respectively for the gate voltages
Vg1 and
Vg2, both of which have the same magnitude as reported in silicon nanowires.
[19] As expected from Eq. (
6), we obtain the crossover temperatures of
TA1 = 203 K and
TA2 = 202 K respectively for the gate voltages
Vg1 and
Vg2, which are consistent with the experimental observation in the conductance platform A of Fig.
2(d). The result shows that the crossover temperature
TA from NNH to VRH is independent of the gate electric field, which may result from the stronger interactions of the thermally activated electrons.
In order to understand the transition of the conductance at temperature TC, we provide the Arrhenius plot of conductance in ln G vs. (1/T)1/2 scales by linear fitting in Fig. 4. With the temperature decreasing below 150 K, the activation energy extracted from Fig. 2(d) is further reduced to be about 0.84 meV and 0.68 meV respectively for the gate voltages Vg1 and Vg2, by which the VRH at the vicinity of Fermi level EF is associated with the long-range Coulomb interaction of the dopant-induced QD array. If the energy of each dopant site depends on distribution of other charged dopants, the electrons localized on the states near the Fermi level will be redistributed. The long-range Coulomb interaction energy
could reduce the DOS for the electrons localized on the states of QDs around the Fermi level, in which the number of states is given by
. In order to estimate the DOS near the Fermi level, ES-type VRH model considers that the Coulomb interaction energy
is comparable to the activation energy
. Then, the DOS near the Fermi level is given by[21]
 | |
whose curve is shown in the inset of Fig.
4. At the crossover temperature
TC for the VRH transition from Mott-type law to ES-type law, the activation energy
Ea should be equivalent to the Coulomb interaction energy
, i.e.,
The DOS at the condition of
would remain constant as
g0 in Mott VRH, which is given by
Then, the hopping distance
of ES-type VRH can be estimated by the maximum hopping probability
[33] of Eq. (
3) as
in which the localized length
in ES hopping regime indicates the average modulation length of long-range Coulomb interaction between the dopant-induced QDs in the conduction channel. In order to determine the localization length
, we substitute Eqs. (
8) and (
10) into Eq. (
3) to get the ES-type VRH conductance expression
The temperature factor
TES is related to the localization length
by
The temperature factor
TES extracted from the slopes is 34.8 K for the gate voltage
Vg1 and 31.4 K for the gate voltage
Vg2. Therefore, the localization length
for ES-type VRH can be obtained from Eq. (
12) to be about 323 nm and 358 nm respectively for the gate voltages
Vg1 and
Vg2. The result indicates that the localization length
is enhanced with the increase of the gate voltage due to screening of the trapping potentials of the dopant atoms. Considering the expressions from Eq. (
7) to Eq. (
10), we finally obtain the expression of
TC,
which is determined by the constant DOS
g0 and the localized length
aES. As a result, the crossover temperature of the VRH conductance transition from Mott-type law to ES-type law can be estimated to be
TC1 = 85 K at the gate voltage
Vg1 and
TC2 = 126 K at the gate voltage
Vg2. Both the two crossover temperatures
TC1 and
TC2 are consistent with the experimental observation of conductance platform B of Fig.
2(d). The crossover temperature
TC2 of 126 K is much higher than
TC1 of 85 K, which may result from the stronger Coulomb interaction at the gate voltage
Vg2 of the
gm valley.